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Abstract 

Random fluctuations of the shot-noise power in disordered graphene nanoribbons are studied. In particu¬ 
lar, we calculate the distribution of the shot noise of nanoribbons with zigzag and armchair edge terminations. 
We show that the shot noise statistics is different for each type of these two graphene structures, which is a 
consequence of presence of different electron localizations: while in zigzag nanoribbons electronic edge states 
are Anderson localized, in armchair nanoribbons edge states are absent, but electrons are anomalously lo¬ 
calized. Our analytical results are verified by tight binding numerical simulations with random hopping 
elements, i.e., off diagonal disorder, which preserves the symmetry of the graphene sublattices. 
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1. Introduction 

It is widely recognized that time dependent cur¬ 
rent fluctuations due to the discreteness of the 
electrical charges-shot noise power-provide further 
physical information of an electronic system than 
other transport quantities such as the conductance. 
For instance, the shot noise takes into account the 
Pauli principle and it can reveal electron correla¬ 
tions. Markus Biittiker and collaborators recog¬ 
nized the importance of the shot noise for the un¬ 
derstanding of the problem of quantum electron 
transport and made major contributions to this 
topic. 

The shot-noise power can be studied by using a 
scattering approach to quantum transport. Within 
this framework, Biittiker found that the shot-noise 
power spectrum P, in the zero frequency limit and 
zero temperature, can be written as [i|, 0 

N 

P = 2eVGo^r„(l-T„), (1) 

n—1 
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where Go(= 2e'^jh) is the conductance quantum, 
V the applied voltage, while the T„’s are the trans¬ 
mission eigenvalues of the Hermitean matrix iff, t 
being the N x N transmission matrix. If there were 
no correlations among electrons, the shot noise is 
given by the Poisson value Pp = 2eVGo'^^Tn. 
The shot noise power has been extensively studied 
both experimentally and theoretically in small elec¬ 
tronic devices such as quantum wires and quantum 
dots, in which quantum coherence is preserved. The 
literature on this topic is very extensive, we thus re¬ 
fer the reader to the review articles Refs. ii. 

In general, electron correlations reduce the shot 
noise respect to the case of fully uncorrelated elec¬ 
trons. The Fano factor F measures that suppres¬ 
sion of the shot noise and it is defined by the ratio 
F = (P)/(Pp), where the brackets indicate energy 
or ensemble average. Using random matrix theory 
to quantum transport, it has been predicted that 
the Fano factor takes the value 1/3 for disordered 
quantum wires in the diffusive regime limit [^, Q , 
while for ballistic chaotic quantum dots F = 1/4, in 
the limit of large number of channels supported by 
the leads attached to the dots 0, [1|. Thus, univer¬ 
sal values of the shot noise suppression have been 
predicted for both transport regimes. 

With respect to graphene, several electronic 
properties have been intensively studied since its 
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Figure 1: Disordered graphene nanoribbons (shaded areas) 
of length L with zigzag (a) and armchair (b) edges with per¬ 
fect graphene leads (non-shaded areas) attached. The width 
W of the nanoribbons has been fixed: W = lla./\/3 and 
5a for zigzag and armchair nanoribbons (a being the lattice 
constant a 2.46A), respectively. With those widths, the 
attached perfect leads have gapless metallic band structures. 


discovery in 2004 and the shot noise is not an ex¬ 
ception (for a review, see and [l(il|V For pris¬ 
tine graphene structures whose lengths are shorter 
than their widths, it has been predicted a 1/3 shot- 
noise suppression, at the Dirac point 11| . This sup¬ 
pression value coincides with the Fano factor value 
for disordered normal-metal wires in the diffusive 
regime. 

Although pristine graphene structures show 
many interesting electronic properties, in real 
graphene-based devices those pro pert ies may be af¬ 
fected by t he p resence of disorder 13, 12. 13. 14 M 
nil . [T^ . [isl . [13 . . In general, different sources of 

disorder can be present in graphene such as ripples, 
vacancies, adatoms, or distortions of the lattice pro¬ 
duced by interactions with the substrate [U 22 1, 
however, even sus pen ded graphene structures are 
not free of defects Hi ill. 

In particular, effects of the presence of disor¬ 
der on the shot noise in graphene sheets have 
been experimentally and theoretically investigated 


li 25 . 26, 27, li It has been found that the 


Fano factor is affected by the strength of the dis¬ 
order, as well as the length-to-width ratio of the 
graphene sample. Most of those studies, however, 
have been concentrated on wide geometries and 
models of disorder that break the chiral symmetry 
of the graphene sublattices. 

Here, we are interested in the properties of the 
shot noise at the Dirac point in finite disordered 
graphene samples, i.e., disordered nanoribbons, in 
which the edge terminations as well as the symme¬ 
try of the graphene sublattices play a crucial role in 
the electronic properties. We thus investigate the 


random fluctuations of the shot noise power for the 
two different edge terminations: zigzag and arm¬ 
chair, Fig. [TJ For both terminations, we consider 
the so-called off-diagonal disorder (random hopping 
connecting the two graphene sublattices) in order to 
preserve the chiral symmetry of the graphene sub¬ 
lattices. This kind of short-range disorder might 
model distortions like ripples in the graphene lat¬ 
tice. Actually, experimental evidence of short range 
disorder in graphene has been observed [H |3l|, HI . 
Therefore, here we are interested in studying the ef¬ 
fects of the nanoribbon edge terminations and the 
presence of disorder on the statistical properties of 
the shot noise power. 

The most interesting properties of graphene are 
found at low energies, i.e., near the Dirac point 
where a linear dispersion relation holds and an 
analogy with relativistic massless particles has at¬ 
tracted much attention. Also, the lattice symme¬ 
try of graphene or chiral symmetry, resulting in 
a symmetric energy spectrum around the Fermi 
energy plays an important role in the description 
of the electronic properties; hence, we are inter¬ 
ested in calculating statistical properties of the shot 
noise near the Dirac point. We recall that in clean 
graphene nanoribbons the band structure is deter¬ 
mined by the edge termination and the width of the 
nanoribbons. Here we consider disordered nanorib¬ 
bons with attached perfect leads whose band struc¬ 
ture is metallic. In the next section we introduce 
a statistical model to describe the random fluctua¬ 
tions of the shot noise power, Eq. O- 


2. Statistical Model 


The statistical properties of the shot noise power 
in disordered graphene nanoribbons can be ana¬ 
lyzed through the statistical properties of the trans¬ 
mission eigenvalues T„, as we can see from Eq. (ED- 
A well established theoretical framework to study 
the transmission statistics of disordered systems is 
the one-parameter scaling approach to localization 
and random matrix theory 341. Our statis¬ 
tical analysis of the shot noise is based on that 
scaling approach to Anderson localization and a re¬ 
cent extension to the case of anomalous localization 
( 35 I HI • In general, the manner in which electrons 
are localized determines the statistical properties of 
quantum transport. Therefore, the statistics of the 
transmission, or conductance, for standard (Ander¬ 
son) and anomalous localizations is different. Eor 
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instance, the average of the conductance of disor¬ 
dered systems in the presence of Anderson localiza¬ 
tion decays exponentially with the system length, 
whereas the conductance average has a power-law 
decay with the length for disorder systems with 
anomalous localization. Both kind of localizations 
have been investigated in disordered zigzag and 
armchair nanoribbons [s^. 

On the one hand, near the Dirac point, only a 
single transmission eigenvalue contributes to the 
electronic transport. Therefore, since we are in¬ 
terested in the shot noise near the Dirac point, we 
only need to consider one transmission eigenvalue 
and the shot-noise power in Eq. CD is reduced to 
P = 2eVGoT{l — T). For convenience, we intro¬ 
duced the dimensionless shot-noise power S defined 
as: 


Thus, along this work we will be interested in de¬ 
scribing the statistical properties of the shot noise 
power S, as given by Eq. CD- 

On the other hand, within the scaling approach 
to localization and random matrix theory, the dis¬ 
tribution of the transmission P{T) is given by 
the solution of the so-called Mel’nikov’s equation, 
which is an evolution equation of P(T) with the 
length of the disordered system ■ The 

exact solution of the Mel’nikov’s equation is known 
in terms of quadratures. To simplify the calcula¬ 
tions and provide analytical expressions for shot 
noise distribution, we find convenient to consider 
an approximation to the exact solution ^ : 


PsiT) = 


' acosh(l/Vr) „_,-iacosh’^a/VT) 


,(3) 


where C is a normalization constant and s is the ra¬ 
tio L/l, L being the length of the system and I, the 
mean free path. The above distribution, Eq. CD, 
has been verified in a number of numerical simula¬ 
tions for any practical value of the disorder strength 
(measured by the value of s) as well as in microwave 
experiments [d^ . Additionally, the value of the pa¬ 
rameter s can be obtained from the numerical or 
experimental data through the linear dependence of 
the average (— InT) with L: {— InT) = L/l{= s). 

We point out that the distribution given in Eq. 
CD is appropriate for systems in which the pres¬ 
ence of disorder leads to an exponential localization 
of electron wavefunctions (Anderson localization) 


with the distance r: lipl ~ 6 “''''’, 7 being a con¬ 
stant. This exponential decay has been experimen¬ 
tally and theoretically studied in several different 
disordered systems [dSj. 

The presence of disorder, however, can lead to a 
different electron localization, or anomalous local¬ 
ization, in relation to the above-mentioned stan¬ 
dard exponential decay. For example, in one¬ 
dimensional disordered systems, at the center of 
the band, the wavefunction decays as \ip\ ^ e“^’' 


with a = 1/2 [44 l45l| . For disordered armchair 


nanoribbons, it has been also shown that electrons 
are anomalously localized (a = 0.69) [s^- 

Recently, it has been proposed a generalization 
of the single scaling approach to describe the trans¬ 
mission through disordered systems with anoma¬ 
lous localization [s^. In this case, the distribu¬ 
tion of the transmission is given in terms of long¬ 
tailed probability density functions or Levy-type 
distributions qa{x), where a is the power-tail ex¬ 
ponent of the density function, i.e., for large x, 
Qaix) ^ l/x^"*"". The model is summarized by the 
following equation which gives the distribution of 
the transmission in terms of two quantities: the av¬ 
erage (InT) and a: 


pOO 

PiAT)= / PsiT)q^{z)dz, 
Jo 


(4) 


where we have defined ^ = (—InT). Ps{a,i,z)iT) 
is given by Eq. with s replaced by the 

function s{a,^,z) = ^/( 2 z“/q), where la = 

1/2 z~°‘qadz. The density function qa{z) is 
part of a family of probability densities commonly 
known as a-stable distributions. We remark that 
only two quantities (^ and a) determine the dis¬ 
tribution of the transmission. Also, we point out 
that Eq. Q predicts a nonlinear behavior of the 
logarithm of the transmission with the length L: 
(InT) oc T“, which is in contrast to the linear be¬ 
havior expected for the standard Anderson localiza¬ 
tion. The model summarized in Eq. has been 
applied recently to describe the transmission of mi¬ 
crowaves in disordered waveguides |^| 

In the following two sections, we will apply Eqs. 
0 and 0 to calculate the distribution of the shot 
noise power S, Eq. 0 , in zigzag and armchair 
nanoribbons, respectively. The theoretical predic¬ 
tions will be compared with numerical simulations 
using a standard tight-binding Hamiltonian model: 


^ +c]ci), 

<i,j> 


(5) 
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where i and j are nearest neighbors and c| (c^) 
is the creation (annihilation) operator for spin¬ 
less fermions, while the hopping elements tij 
between the two graphene sublattices are ran¬ 
domly obtained from the distribution p{t) = 1/wt 
with exp(—r(;/2) < t < exp(w/2), w being the 
strength of the disorder which is fixed to 1. This 
short-range disorder models random distortions in 
graphene sheets without breaking the symmetry of 
the graphene lattice. The numerical simulations 
were performed near the Dirac point: E = 10“® (in 
units of the hopping energy of the perfect leads) for 
zigzag and armchair nanoribbons of widths W = 
11 a/ a/s and W = 5a (a being the lattice constant 
a ~ 2.46A), respectively. The transmission eigen¬ 
values in Eq. m were calculated by using a recur¬ 
sive Green’s function and, for all cases shown in this 
work, the shot noise statistics were collected from 
an ensemble of 2 x 10"'’ disorder realizations. 


single parameter s, which can be extracted from ex¬ 
perimental or numerical data, as we have pointed 
out. 

In order to verify the above result, we com¬ 
pare the theoretical distribution, Eq. with 

the tight binding numerical simulations, described 
previously. In Fig. [51 the numerical distributions 
(histograms) and Ps{S) (solid line) are compared 
for two different strength of disorder, measured by 
the parameter s. The value of s is obtained from 
the numerical data and it is plugged into Eq. ®. 
Figs. Ha) and Ht>) correspond to zigzag nanorib¬ 
bons which are characterized by the average dimen¬ 
sionless conductance (G) = 0.46 and (G) = 0.04, 
respectively. The Fano factor values for those cases 
are F = 0.37 for Fig. Ha), while F = 0.6 for the 
case shown in Fig. H^)- We can observe a good 
agreement between the model [Eq. (H] and the nu¬ 
merical simulations. 


3. zigzag nanoribbons: shot noise via edge 
states 


Firstly, we consider the case of disordered zigzag 
nanoribbons, Fig. Ha)- In pristine zigzag nanorib¬ 
bons, it is well known the existence of edge states, 
which are perfectly transmitted. In the presence 
of disorder, those states remain at the border of 
the nanoribbons, but they are (Anderson) localized 
[^ . On the other hand, since near the Fermi en¬ 
ergy only a single channel contributes to the trans¬ 
mission, we can use the distribution function in 
Eq. (H to calculate the distribution of the shot- 
noise power given by P{S) = {S{S — T(I — T))). 
Thus, using Eqs. m and it is straightforward 
but lengthy to perform the average indicated with 
brackets. The final expression for the distribution 
of the shot noise power is 


Ps{S) 


c 


V / ^aCOSh ^acosh^ (^/\/‘^+) 


/aCOSh _g ^acosh^ (l/y^lsT) 


( 6 ) 


where G is a normalization constant and S± = 

(1 ± a/ 1 + 4S')/2. We notice that the statistical 
properties of the shot noise are determined by the 
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4. Armchair nanoribbons: shot noise via 
anomalously localized states 

We now consider the case of nanoribbons with 
armchair terminations. Fig. Hb)- have men¬ 

tioned, evidence of the presence of anomalous lo¬ 
calization has been found when the hopping ele¬ 
ments of the tight binding Hamiltonian [Eq. (H] 
have random fluctuations, which may model ran¬ 
dom distortions in the graphene lattice. As in the 
case of zigzag nanoribbons, near the Fermi energy, 
only a single channel contributes to the electronic 
transport and the shot noise distribution can be ob¬ 
tained by performing the average {5{S —T{l — T))). 
In the present case, however, the average indicated 
with brackets is performed with the transmission 
distribution given by Eq. 0. The a parameter of 
qa{z) in Eq. (H that characterize the anomalous lo¬ 
calization of our disordered armchair nanoribbons 
has been numerically obtained (a = 0.69) in (s^ . 
This value of a does not depend on the length, 
width, and strength of the disorder, according to 
the numerical simulations. On the other hand, un¬ 
fortunately, there is no a close analytical expression 
for the Levy-type distribution qa{z) with a = 0.69; 
thus, we express the shot noise distribution for arm¬ 
chair nanoribbons in terms of quadratures: 

P^AS) = J [Ps{S+) + Ps{S.)] qc,{z)dz, (7) 

where Ps{S) is given by Eq. (H with s replaced 
by s = A{‘^z°'Ia)- As we have defined previously: 
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Figure 2: Distribution of the shot noise power for disordered 
zigzag nanoribbons. Panels (a) and (b) show the numeri¬ 
cal (histograms) and theoretical distribution (solid line) for 
disordered nanoribbons with dimensionless conductance av¬ 
erage {G) = 0.46 (s = 1.03) and {G) = 0.043 (s = 6.21), 
respectively. The value of the Fano factor is 0.37 for panel 
(a), while F = 0.6 for panel (b). A good agreement is seen 
between theory and numerical simulations. 


^ = (-InT), la = 1/2/“z “fcd-z, and Si = (1± 

vms)/2. 

We now compare our expression in Eq. © 
with numerical simulations of disordered armchair 
nanoribbons. In Fig. [3]we show the numerical (his¬ 
togram) and theoretical (solid line) shot noise dis¬ 
tributions for two different values of C The values 
of f are extracted from the numerical data and they 
are plugged into Eq. ©. We may also character¬ 
ized the nanoribbons by the average dimensionless 
conductance: in Fig. |3Ka), (G) = 0.46, while in Fig. 
|3Ib) (G) = 0.15. The values of the Fano factor of 
these two cases are F = 0.27 for panel (a), while 
F = 0.35 for panel (b). We can observe that our 
model describe correctly the trend of the shot noise 
distributions. As in the previous case of disordered 
zigzag nanoribbons, we point out that no parame¬ 
ters have been adjusted in our theoretical results. 

Finally, we would like to contrast the statistics 
of the shot noise power for both type nanoribbon 
terminations; hence, in Fig. 2] we show the shot 
noise distributions for both armchair (dashed line) 
and zigzag (solid line) nanoribbons. For both cases, 
we have chosen nanoribbons with the same aver¬ 
age of the dimensionless conductance ((G) = 0.46), 
shown in Figs. UKa) andEKa). For this value of (G), 
the Fano factor values are: F — 0.37 for zigzag 
nanoribbons, while F = 0.27 for armchair nanorib¬ 
bons. From Fig. |4l the main differences in the dis¬ 
tributions are seen at small values of the shot noise. 
We can also observe that the probability of having 
small values of the shot noise is larger in armchair 
nanoribbons than in zigzag nanoribbons. This can 
be explained by the larger fluctuations in the trans¬ 
mission eigenvalues in armchair nanoribbons, as a 
consequence of the presence of anomalous localiza¬ 
tion, which lead to large fluctuations of the shot 
noise power. Thus, Fig. H] shows an example of the 
effects of the edge termination on the statistics of 
the shot noise power. 

5. Summary and Conclusions 

We have investigated the statistics of the random 
fluctuations of the shot noise power in disordered 
graphene nanoribbons with zigzag and armchair 
edge terminations, near the Dirac point. Within a 
scattering approach developed by Biittiker and col¬ 
laborators, the shot noise can be written in terms 
of the transmission eigenvalues. We thus apply a 


5 


















P(S) P(S) 


30 




Figure 3: Shot noise distribution for disordered armchair 
nanoribbons. Panels (a) and (b) show the numerical (his¬ 
tograms) and theoretical distribution (solid line) for disor¬ 
dered nanoribbons with transmission average (G) = 0.46 
({—InG) = 1.37) and (G) = 0.15 ({—InG) = 6.37) for pan¬ 
els (a) and (b), respectively. The value of the Fano factor 
is 0.27 for panel (a), while F = 0.35 for panel (b). We can 
observe that the theoretical predictions reproduce the trend 
of the numerical distributions. 



Figure 4: Shot noise distributions for zigzag (solid line) and 
armchair (dashed line) nanoribbons with the same average 
dimensionless conductance {{G) = 0.46). Stronger trans¬ 
mission fluctuations due to the presence of anomalous lo¬ 
calization in armchair nanoribbons increase the probability 
of smaller values of the shot noise power, in relation to the 
zigzag nanoribbons. 


random matrix theory of quantum transport to cal- 
culated the complete distribution of the shot noise 
power for both type of nanoribbon edge termina¬ 
tions. 

Our theoretical model predicts different statisti¬ 
cal properties of the shot noise for zigzag and arm¬ 
chair nanoribbons. Those differences come from the 
fact that electronic edge states in zigzag nanorib¬ 
bons are exponentially localized in space, i.e., elec¬ 
trons are (Anderson) localized, while in armchair 
nanoribbons, electrons are anomalously localized 
or delocalized, in relation to the case of Ander¬ 
son localization. Anomalous localization produces 
stronger random fluctuations of the transmission 
which lead to a different statistics of the shot noise. 
We point out that while for zigzag nanoribbons the 
shot noise distribution depends on a single parame¬ 
ter, for armchair nanoribbons the shot noise distri¬ 
bution is determined by two parameters. Those pa¬ 
rameters are not free in our model in the sense that 
their values are obtained from the numerical simula¬ 
tions, or experimental data, and they are used as an 
input in the analytical expressions. Our theoretical 
predictions have been verified by numerical simu¬ 
lations of graphene nanoribbons with off-diagonal 
disorder, which preserve the chiral symmetry of the 
graphene sublattices. We have found that the value 
of the Fano factor is not universal but it depends 
on the strength of the disorder, which is in contrast 
to the known universal value 1/3 for short and wide 
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clean graphene structures. 

Therefore, to conclude, we have shown that 
nanoribbon edge terminations as well as the pres¬ 
ence of disorder play are relevant role in the proper¬ 
ties of the shot noise power in disordered graphene 
nanoribbons. 
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